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We give an analytical formula for the vacuum polarization of a massless minimally coupled scalar 
field at the horizon of a rotating black hole with subtracted geometry. This is the first example of 
an exact, analytical result for a four-dimensional rotating black hole. 


Quantum field theory in curved spacetime can be used 
to understand a lot of interesting features of black holes 
in a semiclassical approximation, most notably particle 
production near the black hole horizon [1]. The calcula¬ 
tion of vacuum polarization or (for a scalar field) is 
the simplest standard probe of quantum fluctuations in 
a black hole background, and can also be used to under¬ 
stand the symmetry breaking and Casimir effects near 
a black hole. Computation of is also a preliminary 
step in evaluating the stress energy tensor which 

contributes to the backreaction through the semiclassical 
Einstein equation. 

Candelas studied the vacuum polarization of a scalar 
field in the Schwarzschild black hole [2] and was able 
to obtain an analytical expression for {fp') at the hori¬ 
zon. Candelas’ methods extend easily to charged static 
black holes; there have also been numerical studies of vac¬ 
uum polarization of scalar fields on general static black 
hole backgrounds beyond the event horizon (e.g. [3] for 
asymptotically flat solutions and [4] for the asymptot¬ 
ically anti-de Sitter case), and analytical computations 
at the horizon of a black hole threaded with a cosmic 
string [^. The case of rotating black holes is much more 
challenging. Frolov [6] was able to calculate the analyt¬ 
ical expression for {(jp) only at the pole {9 = 0) of the 
event horizon, and Ottewill and Duffy [7] have provided 
a numerical evaluation throughout the black hole hori¬ 
zon. However so far no one has been able to give an 
analytical formula for {fp) throughout the horizon of a 
four-dimensional rotating black hole. (An analytic ap¬ 
proximation good for fields with large mass is available, 
however [5], and exact results are obtainable in d = 3 
with AdS asymptotics [TO].) 

In this case we will be studying a particular example 
of rotating black holes that exist in ’’subtracted geome¬ 
try” pTHl4] . Subtracted geometry black holes are non 
extremal solutions of the bosonic sector of N=2 STU 


supergravity coupled to three vector multiplets. These 
black holes are obtained by subtracting some terms in 
the “warp factor” of the original black hole metric in such 
a way that the wave equation for a massless minimally 
coupled scalar field becomes separable and analytical so¬ 
lutions are obtainable. This subtracted black hole metric 
effectively places the black hole in an asymptotically con¬ 
ical box and mimics the “hidden conformal symmetry” 
m of the wave equation on rotating black holes in the 
near-horizon, near-extremal, and/or low energy regimes, 
which is a key motivator for the Kerr/CFT conjecture 
(see e.g. m)- The energy density of the matter fields 
in this new geometry falls off as second power of radial 
distance, thus confining thermal radiation. The classi¬ 
cal near horizon properties of the subtracted black hole 
are the same as the original black hole ones; in particu¬ 
lar, the classical thermodynamics of the subtracted black 
hole is analogous to the standard one although loop 
corrections to the horizon entropy differ |18j). 

The horizon vacuum polarization in the static sub¬ 
tracted metric was studied in m- In this letter we shalll 
consider the subtracted geometry of the uncharged rotat¬ 
ing Kerr black hole. We shall see that the special features 
of the subtracted rotating metric, in particular the well- 
defined nature of the thermal vacuum and the solvability 
of the wave equation, allow us to obtain analytical results 
that are unavailable for the standard Kerr black hole. 

The subtracted Kerr metric is given by: 


ds^ =—A ^PG {dt 4-Ad(fy 


4-A^P 4-d9'^ 4-^siiP 6d(p'^'^ . (1) 


with 


X = r'^ - 2Mr + , G = - 2Mr 4- cos^ 9 

2Mar sin^ 9 


A = 


G 


A = 8M^r - cos^ 9 . (2) 
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(The only difference between this metric and the stan¬ 
dard Kerr metric is the form of the “warp factor” A. For 
the explicit form of gauge potentials and axio-dilatons of 
the STU model, supporting this geometry, see [13].) The 
horizons and their surface gravities and angular velocities 
are given by: 


r±= M± v/M 2 - a2 , 


1 


K± = 


2M 
n± = K± 


M 


VM2 - a2 
a 


± 1 


VM 2 - a 2 


( 3 ) 


We switch to co-rotating coordinates with 

the new angular variable being defined by: 


ip = If — n_i_t. 


( 4 ) 


These are adapted to observers co-rotating with the black 
hole at the horizon. A noteworthy feature of subtracted 
geometry is that outside the horizon there is a globally 
defined timelike Killing vector, written as dt in the co¬ 
rotating coordinates This guarantees that there 

are no superradiant modes and ensures the existence of 
a Hartle-Hawking-like vacuum state adapted to the co¬ 
rotating observers. This is different from the case of or¬ 
dinary Kerr black hole, where there is no such Killing 
vector | 22 l |23| and a physical co-rotating vacuum re¬ 
quires enclosing the black hole in a reflective box jT] | 21 | . 
The subtracted Kerr resembles more in this respect the 
Kerr/AdS black hole (Sj- 

The general algorithm we follow for computing the 
horizon vacuum polarization in the Hartle-Hawking 
state starts by defining the Euclidean Green’s function 
Gh{x,x') (in a state regular at the horizon and infinity, 
and where the modes are adapted to co-rotating coordi¬ 
nates). Then we will evaluate —iGh with radial point 
splitting, perform the mode sum, and subtract the co¬ 
variant divergent counterterms. 

After writing the metric in coordinates {t, r, 9, p) we 
perform the Wick rotation setting t = —ir. The metric 
becomes: 

[AAp - i{l + An+)dTf 

. ( 5 ) 


-I- -I- d9'^ ^ sin^ 9 {dp — ifl+dr) 

Y Cr 


On writing the massless minimally coupled wave 
equation and proposing a solution of the form 

obtain straightfor¬ 
wardly a radial equation which, in the re-scaled variable 
X = (r — ^{r+ + r_))/(r_|_ — r_), reads: 


- -"l 

-dx \ 4/ 

l^mn 


d V? 

dx 4{x- 4 ) 


l{l + l) 


Xlmn{x) — 0 , 


( 6 ) 


where 


{dmn 


2 Mn^r_ — a^{AnA + lA) — 4Aamnr_ 


( 7 ) 


Two independent solutions of the equation, respec¬ 
tively regular at the horizon and at infinity, are: 


/ 1 . 2 ) ^ l) ' 

i^ + h) 

where 


X'lmn ~ / (aimnAlmmCin ''t ; (8) 


4n =n + l, =21 + 2, = 


( 1 ) = ^ _ 2 ( 2 ) ^ 4 

cc -I- 4 ’ a: -I- 4 ’ 


{o-lmn, blmn) — ? + 1 + ± 


( 9 ) 


and the symmetry of the hypergeometric function makes 
irrelevant which branch of the square root is chosen. 

The full Green’s function is expanded as 


Gni-ix, x,9,p ; -ir', x', 9', p') = 


in 

271-ro 


E 


^inK{T—T'') 


n——oo 


Y, E Yr{9,p)Yr{9',p')Gr,,Ux,x'), 


( 10 ) 


l—O m— — l 


where tq = r_|_ — r_ = 2y/ M 2 — k = as defined in 
and 

G . (t = 7’(Omin)r(&min) (1) / x (2) / x 

mln\‘^-> ) P^2/ “h 2)r (1 “h ^/ ' 

( 11 ) 

To evaluate the vacuum polarization at the horizon we 
set X = 1/2, a:' = e + 5 (note that this is a dimension¬ 
less regulator e = {r' — r)/ro) and join the points in the 
other directions, calling the resulting Green’s function 
Gh{^,9). All the terms in the sum vanish except n = 0, 
so we are reduced to: 


-iGh{+9) = 




87 r 2 ro ^ ^ (I+ m)l 

l—O m— — l 


r(Z -I- 1 -I- iam)T{l -|- 1 — iam) ^ 


r{2l + l) 

X A ( ^ -|- 1 -|- iam, I + 1 — iam, 21 + 2, 


1 -he 


( 12 ) 


where the parameter a = a/r_|_ takes values between 0 
and 1. We replace the hypergeometric by an integral 
expression using formula 9.111 of [3S], leading to: 


00 L /1 \ I 

-*G^(e,0) = ^5](2Z + l) 5] )vE^[Pz"*(cos0)]^ 


dt 


t{i-t)y 1 


\ C — t J l-t -6 — t 


COS (main A) , (13) 































3 


where A= 

The addition theorem for the associated Legendre 
polynomials is used to compute the sum over m, and 
formula IIL4 from |26j subsequently yields the sum over 
1. This leads, after a change of variables to x = 1 — t, to 
the integral expression 


-iGnie^e) 


K 

To 



dxfe{x ); 


(14) 


explicitly the added coutnerterms. This leads to: 

1 


dw 


sin^ (a sinh 

1 v"^a^{\ + a^) 


+ 


w'^{l + a^v‘^)'^G 2w{l-\-w){l-\- 

1 w^q:^( 1 + a^) Ine 


2e^l^(\a^v'^YG 12(1 + 


(19) 


fe{x) = - 


+2tx+(2-x)x^ 

(x2+e)3 


l+^^^^(yi^^sin"0sin2(f InA) 


3/2 


(15) 


where ~ stands for equivalence as e —0. The second 
subintegral is thus reduced to a finite integral involving 
no regulator, that can be evaluated numerically, plus two 
explicit divergent terms. 

In the first subinterval, we can show that: 


with A = X{t{x)). It is easy to see from numerical evalua¬ 
tion that the leading divergences in the integral as e —>■ 0 
match those provided by the standard counterterms m. 


( 10 ) 


where cr is the halved geodesic distance between the 
points and /r is an arbitrary mass scale. It is more dif¬ 
ficult, however, to obtain an explicit expression for the 
finite result of the subtraction. To make progress we per¬ 
form the following sequence of changes of variables: 



( x{l + e) \ 
{{l-x){x + e)) ’ 


w = sinhw. (17) 


This leads to the more tractable expression for the inte¬ 
gral/^ = fgdxfYx): 


G = 


dw- 


VT~ 


'o [e-I-( 1 -I-e)w 2 -I-sin^(asinh ^ w)] 


1 3/2 


( 18 ) 

where v = sind. The intermediate M-integral expression 
is also obtainable directly from dimensional reduction 
from the Euclidean Green’s function in AdS^xS^, using 
the higher-dimensional embedding of subtracted geome¬ 
try described in [T ^0 

To analyze the small e limit and subtract explicitly the 
counterterms, we set aside momentarily the Vl + e pref¬ 
actor and split the integral in two subintervals. If over 
(0, e^/®) and If over (e^/®, -boo). In the second subinter¬ 
val we can set e to zero, at the expense of an error that 
vanishes as e —> 0. Then we can add and subtract terms 
compensating for the leading divergences at the lower 
limit, take e —t 0 safely in the subtraction, and integrate 


^ We thank Finn Larsen for bringing this point to our attention. 


A/e 


/< = 


dw 


0 [e-b ( 1 -b e)w 2 -b sin^(asinh ^ w)] 
die 


3/2 


A/e 


-b (1 -b e)w^ + 


i 2 (q 2 _|_]^N^ 13/2 ’ 

( 20 ) 


which is expressible (formula 3.163.3 of [5^) in terms of 
the incomplete elliptic integrals of first and second kind, 
F( 7 ,fc) and E{j,k). Here 


7 = arcsin 


I c_ -b c+ 
7 ^ V c_ -b ei /3 


k = 


c+ 


C- + c+ 


and c± are the coefficients appearing in the denominator 
of the integrand when it is factored in a form proportional 
to [{cf — w^){cf + We need the expansions of 

the elliptic functions near ( 7 , fc) = (f,l), which have 
been derived in [5B]. In order to obtain all the divergent 
and finite contributions to If, we need F accurately to 
order 1 and E accurately to order e. This in turns require 
obtaining the argument k accurately to order e and 7 to 
order . The result of this expansion is the following 
expression for the divergent and finite pieces of If'. 


= ^ 2ei/3(i + a2u2)3/2 + evTT^ 

+ 6(lTip2)5/2 ^ + 

-b 0:^(1 -b a^)t'^(ln(8(l -b 


( 22 ) 


There is an additional finite contribution coming from 
the prefactor fX -b e to the integral, which yields when 
expanded a 1/2 multiplied by the coefficient of the linear 
divergence of the integral. The complete result is thus 
expressed as: 


A = // 


2\JX -b cYv^ 


(23) 
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with the first to terms given by (22) and (191 respectively. 
We see that the divergences cancel out, leaving only 
linear and logarithmic divergences that will match those 


of counterterms (161, leaving a finite renormalized result. 


This concludes the computation of the explicit diver¬ 
gent and finite portions of the Green’s function’s coin¬ 
cidence limit. The counterterms ( |16[ ) need now to be 
evaluated as a function of e to the order 0(1). The form 
of cr can be computed from the formulas expressing a in 
terms of coordinate separation: 


(j = -I- Aabc^x°‘ 

+ BabcdA>.x°'Ax^Ax‘^Ax‘^ -!-••• 


(24) 

where A, B are obtained from symmetrized derivatives of 
the metric tensor, as described in 

These expressions are valid in a coordinate system in 
which the metric is regular. We use the Kruskal coordi¬ 
nates for the subtracted geometry that have been derived 
in [^, which take the form {U,V,0,ip) with {—UV) oc 
(r — r_|_) near the horizon. Our radial coordinate sepa¬ 
ration is therefore written as Ax°' = (—S,S,0,0) (with 
S oc y^). After computing cr by this procedure (leading 
to an expression of the form cr = -I- /32e^ -I- 0(e^)) it is 


easy to obtain the Ricci counterterm in (16) because to 
the relevant order 0(e) we have = Rrr<x’''cr’^. 

Once all the counterterms are computed by this pro¬ 
cedure, when expressed in terms of the a parameter they 
take the relatively simple form: 


G div — 


1 + 

647r2 M2 


1 


EVT+”a 2 ^ 


-b a^) Ine 
4(1 + a2u2)5/2 


(—1 + q:^(—4 + + {1 + + a^)v‘^ + Sa^u"^)) 

12(1 + 

(25) 

(plus a term of the form i?(r+, ff) In/c^). Then, absorbing 
some i?-proportional terms into the arbitrary constant /c, 
the final result for the vacuum polarization is: 


= R{r+-,Q) In^i^ + 


1 -b o:^ f 1 

647r2 M2 I 12(1 -b Q!2r;2)5/2 


(1 - a"(-4 + aY9 + 9a" + - 3a"u^)) 

— 3a"(l -b a")" ln(l -b a"?;") 

1 


du> 


[ry2-b ^2 sin" (a sinh ^w)]^^ 

1 r;"a"(l-ba") 


y^(l-b a2u2)3/2 2i(;(l-b w)(l-b a2w2)5/2 


where 


R{r+,0) = 


3a"(l -b a")"r:" 

8 M2(l -b a2'(;2)3/2 


(26) 

(27) 


64 TT^ (02J 
0.7- 



sin(0) 


FIG. 1. Vacuum polarization (without R term) at the horizon 
as a function of u = sin 6, for a = a/r-|- = 0 (full line), a = 0.5 
(dashed line), a = 0.75 (dotted line), and a = 1 (dot-dashed 
line). 


The angular profile for the vacuum polarization, ne¬ 
glecting the arbitrary term proportional to i?, is de¬ 
picted in Figure 1. Notice that in the absence of ro¬ 
tation spherical symmetry is recovered, with its value 
= (7687r^M")“^ matching the result obtained 
in [T9] for the subtracted Schwarzschild black hole. In ad¬ 
dition, the result at the pole takes the form (^!>")r+,e=o = 
(7687r^M^)“^(l -ba")(l -b4a" — a^), agreeing with result 
found in [11] using a non-corotating vacuum state (at the 
pole, the distinction is irrelevant). The dot-dashed plot 
corresponds to the extremal case a = M. 

It would be interesting to compare our results with 
numerical computations of the vacuum polarization in 
the standard Kerr metric (with a mirror in place to de¬ 
fine the vacuum). Our calculation holds for the mini¬ 
mally coupled field, and the numerical results in [7] are 
for the conformal case, so a direct comparison is not yet 
available. We expect our calculations to be easily gen¬ 
eralized to the case of fields with higher spins as well as 
to rotating charged black holes, including multi-charged 
solutions [30H32]- We also expect our methods to be ap¬ 
plicable to the computation of the stress-energy tensor, 
which would open the possibility of using the subtracted 
approximation to study analytically the backreaction for 
rotating four-dimensional black holes. 
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